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Abstract 

LHC constraints on the recently proposed SUSY SO(IO) GUT model with top-bottom-tau 
Yukawa unification are investigated. In this model, various phenomenological constraints are 
in concord with the Yukawa unification thanks to the negative sign of /i, D-term splitting in 
the soft scalar masses and non-universal gaugino masses generated by a non-zero F-term in a 
24- dimensional representation of SU(5) C SO(IO). After discussing the impact of the CP-odd 
Higgs boson mass bound on this model, we provide a detailed analysis of the recent direct 
SUSY searches performed by ATLAS and investigate the constraints on this SO(IO) model. At 
95% confidence level, the lower limit on the gluino mass is found to be 675 GeV. Assuming an 
integrated luminosity of 10 fb~^, this bound may be extended to 1.1 TeV if the right-handed 
down squark is lighter than about 1 TeV. 
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1 Introduction 



In the simplest versions of supersymmetric (SUSY) Grand Unified Theories (GUTs) based on 
the SO(IO) gauge symmetry, in which the Minimal Supersymmetric Standard Model (MSSM) 
Higgs doublets, Hu and Ha, reside in a 10 of SO (10), the Yukawa couplings of top, bottom 
and tau unify at the GUT scale [1]. Unlike the gauge coupling unification, Yukawa coupling 
unification is very sensitive to weak scale threshold corrections which depend on the entire 
SUSY spectrum [21 [3l H]. Moreover, top-bottom-tau Yukawa unification is often in conflict 
with radiative electroweak symmetry breaking (REWSB) [3J. Therefore, the condition of the 
top-bottom-tau Yukawa unification at the GUT scale imposes non-trivial constraints on the 
soft SUSY breaking terms. In particular, it has been shown that the Yukawa unification is 
incompatible with the REWSB in the Constrained MSSM (CMSSM) where the soft masses for 
scalars and gauginos are universal at the GUT scale [3]. 

SO (10) models are rather constrained not only because of the requirement of Yukawa unifica- 
tion but also because all the MSSM sfermions of each generation fit into a single 16 of SO (10). 
Thus, assuming the Minimal Flavor Violation [5l |6], the only non-universality in the scalar 
masses consistent with the SO(IO) gauge symmetry is the mass splitting between sfermions and 
Higgses. This has been shown to be insufficient to make the Yukawa unification compatible 
with the REWSB [7]. However, if the SO(IO) gauge symmetry is spontaneously broken, as is 
generally expected, in the effective theory below the GUT scale, a generic D-teim contribution 
splits the scalar masses residing in the same representation of SO(IO) |8]. Most importantly, the 
D-term contribution splits the masses of the MSSM Higgs doublets and may make the Yukawa 
unification consistent with the REWSB |9j. 

The requirement of the top-bottom-tau Yukawa unification prefers negative values of the 
Higgs-mixing parameter, /xj^ However, a model with the negative /i and the universal gaugino 
masses [10] predicts a negative SUSY contribution to the muon anomalous magnetic moment, 
{g — 2)^, and enlarges the observed discrepancy between theory and experiment. That is why 
most of the studies that assume the universal gaugino masses in the last decade have been 
devoted to the case of the positive /i [HI [12]. In such a case the Yukawa unification can be 
realized but there is a price for it. Namely, the scalar masses are pushed to the multi-TeV 
regime so proper REWSB is obtained at the expense of severe fine-tuning. Moreover, these 
models have serious problems with the FCNC constraints [13] because of a very large value of 
|/iy4j|, in particular with a lower limit on BR(i?s — )■ fi~^fi~) which was recently improved by the 
CMS and the LHCb experiments [H]. 

Last but not least, in that class of models the solutions with the most accurate unification of 
the Yukawa couplings at the GUT scale found in the numerical scans performed in Ref. [Ill 1121 
[15] are characterised by the gluino mass of about 350 GeV Q while the recent ATLAS "b-jet" 

"'^We use the sign convention in which the gluino mass parameter M3 and the gluino mass rrig are positive. 

^ The usually dominant SUSY correction to the bottom mass is proportional to /Ltmg/m~ when m^^ > rrig. 
Assuming the Yukawa unification and /i > 0, this correction drives the bottom mass away from its experimental 
value so gluino is required to be much lighter than the heavier sbottom to suppress this correction. Typically 
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search [T7] excluded the gluino masses up to 450 (570) GeV in a model with "ad hoc" splitting 
of the Higgs masses [H] (D-term splitting of the scalar masses [12] )■ For a recent review of 
SO(IO) models with the third family Yukawa coupling unification see Ref. [H]. 

A remedy to the problems of the SO(IO) models discussed above was put forward recently in 
Ref. [19] where a novel model with negative /i was proposed. In this model the gaugino masses 
are assumed to be generated by an F-term which is a non-singlet of SO(IO) transforming as a 24 
of SU(5) C SO(IO). Such an assumption leads to the gaugino masses with the following pattern 
at the GUT scale: Mi : M2 : M3 = — 1 : —3 : 2 [20j. The main phenomenological success of 
the model is its consistency with all experimental constraints including the {g — 2)^ and the 
BR(6 — 57) which are satisfied thanks to the positive sign of the product /XM2 [21]. Another 
advantage of the model is the prediction of a rather light SUSY spectrum which can be tested 
at the LHC In particular, the gluino mass is predicted to be in the range of about 500 — 700 or 
900 — 1550 GeV. It is not unreasonable to expect that some of the parameter space corresponding 
to such a light gluino has already been excluded by the LHC. However, the LHC experiments 
interpret their data mainly in the CMSSM and some simplified models and the obtained mass 
limits for SUSY particles cannot be applied directly to other models. There are several studies 
that reinterpret the LHC SUSY searches in terms of the SUSY models other than the CMSSM. 
The LHC constraints on gauge mediation [22l 123] , anomaly mediation [21], phenomenological 
SUSY models [251 126] , SUSY models with compressed spectra [27] and light third generation 
models [281 ISHl [301 1311 IMl [33] have been studied. The constraints on models with the top- 
bottom-tau Yukawa unification in the context of SUSY SU(4)c x SU(2)i^ x SU(2)ij [34] have 
also been investigated [35]. In the present paper we study the constraints and implications of 
the recent LHC data on the SO(IO) model proposed in Ref. [19]. We focus on three ATLAS 
analyses of the final states with large missing transverse momentum and no leptons: the "0- 
lepton" search with 1.04 fb"^ [36|, the "b-jet" search with 0.83 fb"^ [H] and the "multijets" 
search with 1.34 fb~^ of data [37] . 

The paper is organized as follows. In section [2] we describe in some more detail the SO(IO) 
model proposed in Ref. [19]. In section |3] we summarise the ATLAS searches for squarks and 
gluinos. In section |4] we present our methodology of simulating the ATLAS searches and its 
validation. Interpretation of the ATLAS searches in the SO (10) model proposed in Ref. [19] is 
given in section |5] . We conclude in section [6| In the Appendix we give additional information 
on the validation of our simulation of the ATLAS searches. 

2 Description of the SO (10) model 

We start with a brief review of the SO(IO) model proposed in Ref. [19]. The first crucial 
feature of this model is that the Higgs-mixing parameter is negative. The negative sign 

this implies very light gluino in accord with findings of Ref. [TTJ [TH [H]. However, the necessity of very light 
gluino in the Yukawa-unified models with > was questioned recently in Ref. [TB] where new solutions with 
exact Yukawa unification were found with the gluino mass as large as 1.4 TeV. 
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of /i is phenomenologically preferred because it allows for the negative threshold correction 
to the bottom mass [21 El H] (as required by the Yukawa unification) without the need to 
push up the scalar masses to multi-TeV regime. Even though negative fi is preferred by the 
Yukawa unification, it leads to discrepancy with the experimental results for the [g — 2)^ 
and the BR(6 — )■ sj) if the gaugino masses are universal. As explained in detail in Ref. |19j . 
phenomenologically acceptable models require /iM2 > 0. In order to satisfy this condition, 
the gaugino masses are assumed to be generated by a non-zero F-term in a 24-dimensional 
representation of SU(5) C SO(IO). This leads to the following pattern of the gaugino masses 
at the GUT scale [20]: 

(Mi,M2,M3)= (^-^,-|l^Mi/2. (1) 

It should be stressed that even though non-universalities in the gaugino masses are introduced, 
only one free parameter, M1/2, governs the gaugino sector. 

The soft scalar masses also need to be non-universal to achieve the Yukawa unification and 
light sparticle spectraj^ It is assumed that the D-term of the U(l)v (SO(IO) D SU(5) x U(l)v) 
vector superfield acquires a vacuum expectation value at the GUT scalej^ The scalar masses, 
then, can be written as 

"^L = ^lo + 2^ , 



ml^L=ml,-3D. (2) 

Such a contribution generically arises in the effective theory below the GUT scale as a conse- 
quence of spontaneous breakdown of the U(l)v at the GUT scale [8]. Finally, it is assumed 
that the soft trilinear couplings have a universal value Aq at the GUT scal^ 

It was shown in Ref. [9] that for the universal gaugino masses the top-bottom-tau Yukawa 
unification can be realized only ii D > and mio > w^ie- The same conditions are preferred in 
the SO(IO) model studied in Ref. [19] where the gaugino masses are given in Eq. ([!]). Moreover, 
it was noted in Ref. [19] that the Yukawa unification consistent with the REWSB requires 
specific correlations between parameters D, mio and Aq. Provided that these correlations hold. 



^ In Ref. [38] Yukawa-unified solutions with, gaugino masses given in Eq. ([T]) and universal scalar masses were 
found. However, the price for the universality in the scalar sector is very heavy spectrum of squarks and gluinos 
inaccessible at the LHC. Moreover, such solutions predict too small SUSY contribution to the {g — 2)^ to be 
consistent with the experimental data. In the present paper we demand that the constraint from the (g — 2)^ 
is satisfied at least at 2a level. 

^ In Ref. [39j a model with /i < and non-universal gaugino masses given in Eq. ([l|) was considered in which 
sfermion masses are assumed to be universal and "ad hoc" Higgs mass splitting is introduced which explicitly 
breaks SO (10) gauge symmetry. 

^ This assumption requires the existence of a singlet F-term which dominates the soft trilinear couplings but 
gives a subdominant contribution to the gaugino masses. The phenomenology of models without the singlet 
F-term (and with the non-universal A-terms) will be investigated in a separate article. 
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the top-bottom-tau Yukawa unification can be realized for a very wide range of M1/2 and mig. 
Non-trivial constraints on the model are obtained from the BR(6 — )■ 57), the {g — 2)^ and the 
upper bound on the relic abundance of the lightest neutralinos, ^dm^^- Demanding that the 
predictions of the model for these observables are consistent with the experimental data, the 
assumption of Yukawa unification leads to rather definite prediction for the MSSM spectrum. 
This is mainly due to the tension between the BR(6 — )■ 57) and the [g — 2)^: the former 
prefers heavy SUSY spectra while the latter the light ones. On top of that, the WMAP bound 
|10] on the f^DM^^ can be satisfied only in some corners of the parameter space where various 
mechanisms which enhance the annihilation cross-section of the LSP are at work. It was found 
in Ref. [19] that in the part of the parameter space where the Yukawa unification is consistent 
simultaneously with the BR(6 — )■ S7) and the {g — 2)^, the gluino mass is predicted to have any 
value below about 1550 GeV. However, in the region of the parameter space where the gluino 
mass is below about 500 GeV and between about 700 and 900 GeV, there is no efficient LSP 
annihilation mechanism that could reduce the thermal relic abundance of the neutralinos to 
acceptable values. 

For the gluino mass between about 500 and 700 GeV WMAP bound on the ^dm^^ can be 
satisfied due to resonant annihilations of the LSP through the Z boson or the lightest CP-even 
Higgs exchange. Throughout this paper we will refer to the region of the parameter space with 
the efficient LSP annihilation via Z or as the light gluino region. In light of the recent 
LHC data it is not clear whether the model points with such a light gluino are still allowed. In 
Ref. [3^ ATLAS collaboration presented interpretation of the 0-lepton search in the simplified 
model in which all the sparticles except for the gluino, degenerate squarks of the first and 
second generations and the massless neutralino are decoupled. In such a model gluino with 
mass below 700 GeV has been excluded at 95% C.L. for the squark masses up to 2 TeV, with 
the limit increasing to 1075 GeV for equal mass squarks and gluinos. In the light gluino region 
of the SO(IO) model of Ref. [19], the masses of ur, ul and (which are degenerate to large 
extent) are in the range 1.1-1.6 TeV, while the mass of dn (which is always lighter than the 
other first generation squarks due to the D-term splitting of soft scalar masses) is between 
about 900 and 1300 GeV. Applying limits on the gluino and squark masses in the simplified 
model to that SO(IO) model one would exclude the light gluino region. However, the limits 
on the gluino and squarks masses in the simplified model described above should be treated as 
an useful indication to what extent the gluino and squark masses can be constrained with the 
available LHC data. For realistic models, these limits will be typically weaker. One of the main 
purposes of the present paper is to check whether the light gluino region of the SO (10) model 
of Ref. [19] is still allowed by the recent LHC data. 

There are at least two reasons why the limits from the 0-lepton search on the squark and 
gluino masses in the SO (10) model of [19] may be substantially weaker than in the simplified 
model. First, the winos and Higgsinos are always lighter than the gluino at the EW scalej^ 
Therefore, the gluino usually decays to the chargino or neutralino other than the bino-like LSP, 



6 M3 2M2 12Afi at the EW scale. 
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unlike in the simplified model where the gluino decays to qqXi through an off-shell squark. 
Second, in some part of the parameter space the right-handed sbottom is lighter than the 
gluino. For such model points gluino decays exclusively to bi (associated with a b-jet) which 
subsequently decays most of the times to the Higgsino-like neutralino or chargino unless these 
decays are kinematically inaccessible. Longer gluino decay chains result on average in smaller 

of individual jets and such events are less likely to pass the selection cuts. 

In the second region with the gluino mass between about 900 and 1550 GeV, three mecha- 
nisms enhancing the LSP annihilations may be at work: t-channel exchange of the right-handed 
sbottom, resonant annihilation through the CP-odd Higgs, A^, and coannihilations with the 
stau. The first aforementioned mechanism requires the sbottom to be lighter than about 200 
GeV. The Yukawa-unified solutions with such a light sbottom have been found in the numerical 
scan performed in Ref. [19] but only in the part of the parameter space where the LSP mass is 
below 100 GeV. Such solutions are excluded at 95% C.L. by the DO search at the Tevatron [H]. 
The Yukawa-unified solutions, in which the LSP annihilation via the exchange dominates 
the total annihilation cross-section, were ruled out very recently by the LHC direct A^ — )■ rr 
searches. This follows from the fact that the LSP annihilation via the A^ exchange in the 
SO(IO) model of Ref. [19] can be large enough only when m^o < 400 GeV and tan /3 between 45 
and 50, while the ATLAS [42] and CMS [43] analyses excluded m^o below 400 (450) GeV for 
tan/3 ^ 45 (50) at 95% C.L.. The limits on the (m^o,tan/3) plane do not depend significantly 
on the other MSSM parameters jH], and these limits are also applicable to that SO(IO) model 
in a rather good approximation. 

The regions of the parameter space where the WMAP bound is satisfied due to the stau 
coannihilation are not constrained by any experiment so far. The values of the gluino mass 
found in the stau-coannihilation region are between about 1000 and 1550 GeV, the masses of 
the ul and the dL are in the range 1.3-1.8 TeV, while the ur (d^) are lighter by about 100 (300) 
GeV than the left-handed squarks. In contrast to the light gluino mass region, in the stau- 
coannihilation region, the gluino is always heavier than all the third generation squarks and 
often heavier than the d/j. Most importantly, the bi is always lighter than the gluino by at least 
500 GeV which implies that the sbottom production is comparable or even more significant than 
the gluino and first generation squarks production. Notice that the spectra corresponding to 
the stau-coannihilation region are too heavy to be excluded even if one applies the limits found 
in the simplified model with the gluino, the LSP and only the first generation squarks (based on 
the 0-lepton search) and the simplified model with the gluino, the LSP and the lighter sbottom 
(based on the b-jet search) which typically overestimate the true exclusion limits. Nevertheless, 
we include this region in our analysis to confirm this expectation. 

There is also a part of the parameter space where both the A^ exchange and the stau- 
coannihilation contribute to the LSP annihilation cross-section at similar rate and the WMAP 
bound is satisfied. In this region, the required A^ mass has not been excluded by the current 
LHC searches. However, we should comment that such solutions require ttiao to be smaller than 
about 480 GeV, and those solutions may be excluded when the whole LHC data accumulated 
in 2011 are analysed. The gluino and squark masses corresponding to such solutions are similar 
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Table 1: The cuts used to define the signal regions of the ATLAS 0-lepton analysis ^36j. All the 
dimensionful quantities are given in GeV. The A(/)min is the minimum azimuthal angle between p^^^^ 
and pt of the three leading jets. The number of events observed by ATLAS in each signal region, as 
well as the expected Standard Model background are also displayed. The first uncertainty represents 
the statistical uncertainty, while the second is the systematic one. 



to those for the stau-coannihilation region. 

We should also comment on the fact that the only Yukawa-unified solutions corresponding to 
the gluino mass in the range 900 — 1000 GeV are those with a very light sbottom or rriy^o < 400. 
As we pointed out above, the former are excluded by the Tevatron while the latter are excluded 
by the LHC This is an interesting example of a case in which parameter M1/2 and thus the 
gluino mass can be constrained not because of excessive production of the gluinos in the pp 
collisions but as a consequence of some correlations between the MSSM parameters in a specific 
subclass of MSSM models. 

Another interesting feature of the SO(IO) model is the fact that the lightest CP-even Higgs 
is predicted to be below about 114.5 GeV in the region of the parameter space where the — 
and the b sj are simultaneously satisfied. This is dangerously close to the LEP2 bound |15] 
but taking into account the fact that the uncertainty of this prediction is about 3 GeV jl6] 
this region of the parameter space is still allowed by the experimental data. We should also 
emphasise that the lightest Higgs can be much heavier if the constraint from the {g — 2)^ is 
neglected. In particular, the mass of the Higgs can be about 125 GeV, as suggested by the 
excess observed recently by ATLAS |17] and CMS |18], but the corresponding masses of the 
squarks and the gluino are too heavy to be discovered at the LHC with y/s = 7 TeV. 

3 ATLAS searches for squarks and gluinos 

In order to investigate the LHC constraints on the SO(IO) model we focus in this paper on the 
SUSY searches performed by ATLAS, although one can expect that the limits from the CMS 
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Table 2: The cuts used to define the signal regions of the ATLAS b-jets 0-lepton analysis pJTj. 

searches are comparable. The 0-lepton search performed by ATLAS [36] with 1.04 fb~^ of data 
defines five signal regions with at least two, three or four jets. The main kinematic variables 
used in that analysis are the effective mass, rries, and the magnitude of the missing transverse 
momentum E^^^^ = The is defined as the sum of Eif^^^ and the magnitudes of the 

transverse momentum, = \pt\, of two, three or four leading jets for the signal region with at 
least two, three or four jets, respectively. In the "high mass" signal region, however, the is 
defined as the sum of E^^^^ and of all jets with > 40 GeV. All the signal regions require 
high jets and large E^^^^. The cuts applied in each signal region of the 0-lepton search are 
summarized in Table [H 

In the SO(IO) model of Ref. [H], the right-handed sbottom is the lightest among the squarks. 
In some part of the parameter space, the sbottom is even lighter than the gluino. Therefore, one 
may expect that the sbottoms are significantly produced at the LHC, either directly from the 
proton-proton collisions or indirectly from the gluino decays. In consequence, a lot of b-jets are 
expected in a typical final state. That is why we use also the ATLAS b-jets 0-lepton analysis 
[17] with 0.83 fb~^ of data to set constraints on that SO(IO) model. In that analysis, four 
signal regions have been defined with the requirement of three or more jets and at least one or 
two b-jets, depending on the signal regions. In the ATLAS b-jet analysis the same kinematic 
variables as in the 0-lepton search are used. Since one of the main targets of the b-jet 0-lepton 
analysis is constraining models which predict light sbottoms, the cuts on the mes are weaker 
than in the 0-lepton search. The cuts defining the four signal regions of the b-jet 0-lepton search 
are summarized in Table [2l 

Another analysis which may be sensitive to the SO(IO) model was released most recently 
and focuses on search for new phenomena with large jet multiplicities and missing transverse 
momentum [37] using 1.34 fb~^ of data. Throughout this paper we will refer to that analysis as 
the multijets analysis. In that analysis events with leptons in the final state are discarded. Four 
different signal regions have been defined with the minimum number of jets varying between 
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Table 3: The cuts used to define the signal regions of the ATLAS multijets analysis [37] . 



six and eight. For such a large number of jets, large trigger inefficiencies may occur. In order 
to maintain acceptable trigger efficiency a separation ARjj > 0.6 between all the jets with 
above the threshold for a given signal region is required. The cut on E^^^^ / \/ Ht where Ht is 
the scalar sum of the transverse momenta of all the jets with > 40 GeV and |?7| < 2.8 implies 
that also in that search a large missing transverse energy is required. The cuts applied in each 
signal region of the multijets search are detailed in Table |3} 

The primary target of the multijets analysis is to provide increased sensitivity to models 
with sequential decays to many strongly interacting particles. It therefore fits very well with 
the purpose of constraining the SO(IO) model. The multijets analysis has already been proved 
to be useful in the context of the CMSSM for which constraints in the region where mi/2 <C ttiq 
from multijets analysis are comparable (or even better if mo is between about 1 and 1.5 TeV) 
to the ones from the 0-lepton search [37] . 

For each signal region i, the ATLAS papers provide the information on the observed number 
of events no \ the expected number of the Standard Model background events rS^'' together with 
its error o\^\ These numbers are given in Tables [l, 2,3, The agreement between rio^ and 
can place constraints on SUSY models which predict the expected signal events ni*'' with its 
error alP . In the next section we shall describe our procedure of computing Ti'f' and cxi*^ and 
validate it against the ATLAS exclusion plots for the CMSSM (in the case of the 0-lepton [36] 
and the multijets [37] analyses) and the simplified model with the gluino, the sbottom and the 
lightest neutralino (in the case of the b-jet analysis [TT]). 



4 Methodology and its validation 

Let us now discuss our methodology of simulating the ATLAS searches. We use various publicly 
available codes. First, we compute the SUSY spectra with SOFTSUSY [19] whose outputs 
are passed to SUSYHIT [50] which calculates branching ratios for all SUSY particles. We 
simulate SUSY events using Herwig+-|- [51] where the effects of parton shower, hadronization 
and underling events are taken into account. 10000 events are generated for each model point. 
In the simulation we include only production channels with at least one squark or gluino. We do 
not simulate direct pair productions of the neutralinos, the charginos or the sleptons since they 
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are either negligible or the of jets and the effective mass, rries, corresponding to such events 
are too small to pass the selection cuts. The detector simulation is done using Delphes [52] with 
the ATLAS detector card modified according to the ATLAS analyses, where in particular jets 
are defined using the anti-fcr algorithm with AR = 0.4 and the b-jet tagging efficiency is set to 
50%. Finally, the total SUSY production cross-section, 0"nlO) is calculated at next-to-leading 
order with PROSPINO [53]. 

In order to work out the expected number of SUSY events in the signal region i, ns \ we 
apply the cuts detailed in Tables [l| [2] and |3] to find their efficiency e (i.e. the ratio of the 
number of events that survived the cuts to the total number of generated events). The ni^^ 
is given by ctnlo ^ e*^*-* x x A^^\ where A^^^ is a correction factor which accounts for the 
impact of additional data cleaning due to electronics failure in one of the calorimeters in the 
ATLAS detector during some period of data-taking. We take A = 0.9 for the b-jet analysis and 
A = 0.95 for the 0-lepton and the multijets analyses. 

In order to estimate the exclusion limits we mimic the ATLAS analyses to the extent that 
is allowed by the limited amount of information provided in the ATLAS publications on those 
analyses. In particular, the CLs method [5^ is used which employs a variable CLs defined by: 



CLs = , (3) 

1 -Pb 

to set the exclusion limits on a model at (1 — CLs) x 100% confidence level, which means a 
model is excluded at 95% C.L. if CLs < 0.05. In the above definition ps+b is the p-value of 
the background plus signal hypothesis while is the p-value of the background only hypoth- 
esis. Given the information S*^*-* = {ni^\n^\ai^\al^^) for each signal region i, the p-values are 
computed in the following way [55j : 

(i) 

rXo oo 

= ' = J2 , (4) 

n=0 (i) 

n=no 

where P(n|S*^*)) is the probability of observing n events in the signal region i which is governed 
by Poisson distribution with the expectation value for the observed events: 

A»(S», 5,, 6s) = n«(l + 6,al:^) + n»(l + , (5) 

where the impact of the uncertainties on the expected number of events is accounted for by 
the nuisance parameters 6b and 6s- We have neglected the error on the luminosity since it is 
negligible as compared to the errors we included. Assuming that the nuisance parameters have 
Gaussian probability distribution function, the probability of observing n events is given by: 

/oo poo A*^*' /' \ (i)\n 
d6j d6b- ^e-5(^?+^.^). (6) 

The lower limit of the integration is restricted to keep the number of signal and background 
events independently non-negative. 
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Figure 1: Comparison of our 95% exclusion contours with those of ATLAS in the simphfied 
model in which only the gluino, the lighter sbottom and the LSP with mass of 60 GeV are in 
the low energy spectrum, while all the other sparticles are decoupled. In contrast to the ATLAS 
exclusion contour, our cover also some region of the parameter space where rrig < m^^ + m^. 




The ATLAS papers provide the error on the signal cross-section. It is estimated by vary- 
ing the renormalization and factorization scales in PROSPINO between half and twice their 
default value and by considering the PDF uncertainties provided by CTEQ [56j. The resulting 
uncertainties on the signal cross-section vary somewhat from point to point in parameter space 
but it is reasonable to assume that the error on the cross-section for each production channel 
is constant across the parameter space. In our analysis we assign uncertainties of 10%, 20%, 
30%, 35% for the qq^*\ qg, gg, it* and bb* productions, respectively. 

The ATLAS papers do not provide, however, the details on the other systematic uncertain- 
ties on the signal such as the jet energy scale uncertainty or the b-tagging uncertainty. We 
model these uncertainties by taking a constant value for the signal region i. Our estimation 

of the total systematic uncertainty on the signal is then given by ai*^ = ^ (^'^^oss + ^S*'' • order 
to determine we adjust it manually in such a way that our 95% CL. exclusion contours 
match the official ones presented in the ATLAS publications. 

We find that cXgf'^^^ = 0.6 reproduces the official ATLAS exclusion contours from the b-jet 
analysis [17] in the simplified model containing only the gluino, the sbottom and the lightest 
neutralino with the mass of 60 GeV reasonably welQ as seen from Figure [T| In the official 
ATLAS plot exclusion contours are restricted to the region of the parameter space in which the 

^We do note include in our analysis the 3JA, 3JB and 3JC signal regions of the b-jet search because the 
exclusion plots for the individual signal regions has not been made publicly available so it was not possible to 
validate our simulation of these signal regions. We do not expect that inclusion of these signal regions would 
change our conclusions significantly since the 3JD signal region is the most constraining one in the vast majority 
of the parameter space fTT]. 



11 



gluino mass is larger than the sum of the sbottom and bottom masses. In Figure [T] we present 
extended version of that plot to the region of the parameter space where rrig < m^^ + rrib and 
the gluino decays through the off-shell sbottom into bbxi- 

We also find that ^ = af ^ = = 0.3 and 4"°^^ = a^f'^ = = 0.1 give 

a rather good approximation of the ATLAS exclusion contours for the 0-lepton search, as 
seen from Figure A.l in the Appendix. For the multijets ana lysis we find the best fit for 
f7i55) _ ^(8j55) _ ^(6i80) _ ^(Jjso) _ gj^^ gggj^ from Figure A. 2 in the Appendix. 



5 ATLAS constraints on the SO (10) model 

We are now ready to present resulting constraints on the SO(IO) model from our simulation 
of the ATLAS analyses. We analysed 977 points (out of which 490 are in the light gluino 
region) satisfying the top-bottom-tau Yukawa unification at the level of 10% or better found in 
a numerical scan performed in Ref. [19]. All of these points satisfy the following experimental 
constraints (which were computed with appropriate routines in MicrOmegas |57j): 



12.7 ■ IQ-^^ < Saf^^^ < 44.7 ■ 10"^° {2a) ^ 

2.89 ■ 10"^ < BR(6 ^ S7) < 4.21 ■ 10"^ (2a) P] 

BR(5, ^ fx^fi-) < 1.08 • 10^^ [H] 

^DM/i^ < 0.1288 (3(t) 00] 

m?,o > 111.4 GeV [S] 



Notice that we have included the recent constraint on the BR(i?s — )■ /U+/i~) from the combined 
analysis of CMS and LHCb [H]. However, this observable is much less constraining than the 
lower limit on m^o which we set in our analysis at 400 GeV independently of tan /3 (which 
varies between about 45 and 50 in the region that might be affected by the present constraint 
on niAo). We use slightly relaxed bound on m^o due to its (small) model dependence and the 
theoretical uncertainty in the prediction of rriAo. The mass limits on SUSY particles from LEP 
and the Tevatron including the limit on the sbottom mass discussed in section [2] are also taken 
into account. Due to the fact that the uncertainty in the prediction of the lightest Higgs boson 
mass is about 3 GeV |l6], slightly relaxed LEP2 bound [45] is used|^ Even though we require 
the model points to satisfy only the upper WMAP bound, many of them satisfy also the lower 
WMAP bound on the dark matter relic density. 

The main result of our analysis is that 481 out of 490 points in the light gluino region 
are excluded at 95% C.L.. On the other hand, none of the points are excluded in the stau- 
coannihilation region. This can be seen from Figure [2] where the distribution of the CLs value 

^ Very recently ATLAS 47 and CMS announced their preliminary results of the SM Higgs searches 
with 5 fb~^ of data. Due to downward fluctuation in the data ATLAS extended the lower limit on the SM 
Higgs boson mass to 115.5 GeV at 95% C.L.. We do not take this bound into account because the Higgs boson 
production cross-section in the MSSM is usually smaller than in the SM due to the negative contribution from 
a stop loop to the gluon fusion production rate, while ATLAS excluded only a > 0.9(Tsm for rrih « 115 GeV. 
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Figure 2: Distribution of the CLs value in the mig — M1/2 plane based on the combination of 
the 0-lepton, the b-jet and the multijets searches. For each point in the plot, the signal region 
providing the smallest value of the CLs is chosen. The black crosses correspond to the points 
excluded at 95% C.L.. 
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Table 4: Percentage of the excluded points at 95% C.L. in the Hght gluino region by individual 
signal regions, "comb" in the 0-lepton (multijets) column correspond to a combination of all the signal 
regions within the 0-lepton (multijets) search, while "comb" in the last column corresponds to the 
combination of all three ATLAS searches. 



in the rriiQ — M1/2 plane based on the combination of all searches is presented. For each point 
in this plot, the signal region providing the smallest CLs value is chosen. The model points 
depicted by the crosses are excluded at 95% C.L.. 

In Figures |3] and |4| the distributions of the CLs values based on the individual searches 
and their combination are shown in the m^^ — rrig and m^^ — rrig planes, respectively. The first 
thing to note from these Figures is that the 0-lepton and the multijets analyses provide the 
most stringent constraints on the light gluino region but the b-jet search is only slightly less 
constraining ruling out by itself about 60% of points in this region. In Table |4| fractions of the 
excluded points in the light gluino region by individual signal regions are shown. In the 0-lepton 
analysis, essentially only the "high mass" region is relevant, while in the multijets analysis all 
the signal regions are relevant and complement each other. We should also emphasise the 
complementarity between the 0-lepton and the multijets searches. While individually these 
searches constrain the gluino mass to the values above about 620 and 610 GeV, respectively. 
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Figure 3: Distribution of the CLs value in the m^^ — rrig plane based on: (a) the b-jet search 
only, (b) the 0-lepton search only, (c) the multijets search only, (d) the combination of all 
searches. The black crosses correspond to the points excluded at 95% C.L.. 



their combination excludes all the points with the gluino mass below about 675 GeV. The 
multijets analysis is less sensitive to the part of the parameter space where rrig > m^^ and the 
mass splitting between the gluino and the sbottom is small. In such a case the gluino to sbottom 
decay produces rather soft jets. This implies that passing the cuts with large jet multiplicities 
in such a case becomes less likely]^ Appearance of the soft jets in the final state affects also the 
efficiency of the "high mass" signal region of the 0-lepton search but in a milder way, so some 
of the points (but not all) with small mass splitting between the sbottom and the gluino which 
are not excluded by the multijets analysis are excluded by the 0-lepton search. As a result, all 

^ Note that our setup of the event simulation does not take account of the NLO matrix element correction 
with extra hard jets associated with the sparticle production. Including that contribution may slightly increase 
the efficiency for the events with soft jets from SUSY cascade decays. In this sense, our analysis provides a 
conservative limit. 
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Figure 4: The same as in Figure ISj but in the m^^ — irig plane. 



the points in the hght gluino region which are still allowed are characterized by rather small 
mass splitting between the gluino and the sbottom - smaller than about 100 GeV, as can be 
seen from Figure |3j Another feature of the points not excluded by any ATLAS search is rather 
heavy (Ir. The lower bound on its mass in the light gluino region is found to be about 1.2 TeV, 
as clearly seen from Figure |4j 

It is also interesting to note that all the model points in the light gluino region would be 
excluded if the ATLAS had not observed some small excess over the background (about lo") in 
the high mass signal region of the 0-lepton analysis (see the last column in Table [T]). Therefore, 
if the excess observed by ATLAS is just an upward fluctuation of the background, and not 
manifestation of SUSY events, the light gluino region will be excluded very soon. 

Even though all of the points in the stau-coannihilation region are consistent with the 
present LHC data, one can draw some interesting conclusions from our analysis for this case 
too. First of all, the values of the CLs provided by the b-jet and the 0-lepton search are much 
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Figure 5: Distribution of the expected CLs value based on the 0-lepton search with 10 fb 
of the integrated luminosity. 



smaller than those of the multijets search. This is because in the stau-coannihilation region the 



mass of the lighter sbottom is typically comparable to the Higgsino mas£^° and the 61 — ?■ Xit 
decay channel is closed. The sbottom decaying through this channel may produce as many as 
five jets if decays subsequently to XiW so this channel being closed results in a big loss of 
the cuts efficiency in the multijets analysis. Even if the bi — )■ Xit decay channel is open, as is 
the case for a small number of points in the stau-coannihilation region, the Xi usually decays 
to or tiUt rather than to XiW which results in leptons or at least a smaller number of hard 
jets in the final state making it unlikely that the events would pass the selection cuts used in 
the multijets analysis. Moreover, efficiency of the cuts used in the multijets analysis is reduced 
also by the fact that decays dji usually decays directly to the LSP producing only single jet. 

The 0-lepton and the b-jet searches set complementary constraints on the stau-coannihilation 
region. The b-jet search is the most sensitive (i.e. provides the smallest CLs value) to model 
points with very light sbottom where the total SUSY cross-section is dominated by the sbottom 



production ^ while the 0-lepton search provides the best constraints for models with light dn, 
as seen from Figures |3] and |4j 

We also estimate the expected constraints on the stau-coannihilation region in the SO(IO) 
model with 10 fb~^ of the integrated luminosity. The multijets search will not have enough 



Tliere are two reasons for tliis. One is tliat a requirement of tlie bottom-tau Yulcawa unification positively 
correlates with the gluino mass. The second is that due to large negative contributions to the sbottom mass 
from Z)-term and the rcnormalization group effects (coming from a large positive ^0 required for consistency 
with the 6 — > 57 constraint) the sbottom is light (see a detailed discussion in Ref. |19j). 

In the final stages of preparing this work new ATLAS analysis was released [30] in which the limits on 
the uif^^ — TO^o plane were found that surpass those from DO. In that analysis 2 b-jets are required in the final 
state and the cuts are optimized in order to maximise the sensitivity for models with the light sbottom. That 
ATLAS search excludes the model points in the stau-coannihilation region with the sbottom mass below about 
350 GeV if BR(6i — Xib) = 100% (which is usually the case in this part of the parameter space). 
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sensitivity to exclude any points from this region, as expected from the above discussion. We 
found that the most promising for constraining the stau-coannihilation region is the 0-lepton 
search. It can be seen from Figure [5] that the 0-lepton search can probe the gluino masses in the 
stau-coannihilation region up to 1.1 TeV if the right-handed down squark is lighter than 1 TeV. 
In our estimation of the expected exclusions with 10 fb~^, we rescale the statistical errors with 
the luminosity but keep the relative systematic error constant. This is certainly a conservative 
assumption. Thus the reach of the 0-lepton search will be even better. With this conservative 
assumption, the b-jet search will not be able to probe the stau-coannihilation region. However, 
if the systematic uncertainties on the background were to be reduced significantly, for example 
by data driven-methods of systematic error estimation, and/or the cuts were to be refined 
according to the larger amount of data, then the b-jet search may be complementary to the 
0-lepton search, especially in the part of the parameter space with light sbottom. 



6 Summary and conclusions 

We have investigated the LHC constraints on the recently proposed SUSY SO (10) GUT model 
[T9] predicting the top-bottom-tau Yukawa unification. In this model the non-universal gaugino 
masses arising from SO(IO) non-singlet F-term transforming as 24 of SU(5) C SO(IO), D-teim 
splitting of the scalar masses and the negative sign of fi are assumed. The model is able to 
satisfy all experimental constraints including the {g — 2)^, the b ^ s'-f and predict correct 
thermal relic abundance of the neutralinos. In the region of parameter space of the SO (10) 
model where all phenomenological constraints are fulfilled, relatively light MSSM spectra are 
predicted and there is a good prospect of discovering SUSY particles at the LHC. We found 
that this part of the SO (10) model parameter space tightened considerably due to the recent 
results of the LHC experiments which have not found any sign of SUSY so far, but many 
Yukawa-unified solutions are still consistent with all experimental data. 

The phenomenologically viable parameter space of this model can be classified into three 
regions in terms of the dominant LSP annihilation process: the light gluino region (the dominant 
process is Z boson or light CP-even Higgs exchange), the light region (the dominant process 
is CP-odd Higgs exchange) and the stau-coannihilation region (the dominant process is stau- 
coannihilation.). We have argued that the light ^4° region has already been excluded by the 
recent — )■ rr search at the LHC and we have investigated the impact of the recent direct 
SUSY searches on the light gluino and the stau-coannihilation regions in the SO(IO) model. 

Since the constraints on the gluino and squark masses obtained in the simplified model, 
where only the massless neutralino, the gluino and the first /second generation of squarks are 
assumed in a low energy SUSY spectrum, cannot be directly applied to other models, we 
performed a dedicated analysis to assess the impact of the recent direct SUSY searches on the 
SO(IO) model. We focused on three ATLAS analyses: the 0-lepton search with 1.04 fb~^ |36j . 
the b-jet search with 0.83 fb~^ [TTj and the multijets search with 1.34 fb~^ of data. We 
found the 95% C.L. lower limit of 675 GeV for the gluino mass, substantially weaker than the 
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one obtained in the simplified model. This implies that the light gluino region is excluded 
except for the corner of the parameter space where the gluino mass is between 675 GeV and 
700 GeV, the right-handed down squark mass is above 1.2 TeV and the mass splitting between 
the gluino and the lighter sbottom is smaller than about 100 GeV. We, however, stress that 
the light gluino region would be completely excluded in the absence of the slight excess of the 
events observed in the "high mass" signal region of the ATLAS 0-lepton search. Therefore, this 
region will be excluded very soon unless the observed excess is not just an upward fluctuation 
of the background. 

We also found that the stau-coannihilation region remained unconstrained by the present 
searches. We have estimated the expected exclusion limit on this region assuming the same 
analyses but with 10 fb~^ of the integrated luminosity at the 7 TeV LHC and found that the 
0-lepton search may set the 95% C.L. lower limit on the gluino mass at 1.1 TeV if the right- 
handed down squark is not heavier than about 1 TeV. The b-jet search may also contribute 
to the constraints on the stau-coannihilation region, especially in the part of parameter space 
with light sbottom, if the systematic error in that search is reduced and/or the selection cuts 
are refined. 

In conclusion, the SO (10) model of Ref. [19] is still consistent with all available direct and 
indirect constraints, with good prospects for testing it at the LHC in the near future. This 
should be contrasted with the Yukawa- unified SO (10) models with the universal gaugino masses 
where the most favourable part of the parameter space has already been excluded by the LHC 
experiments. For this reason we believe that the model we have studied in the present paper is 
the best motivated and the most promising SO(IO) GUT model that predicts the top-bottom- 
tau Yukawa unification. Therefore, it is certainly worth further detailed investigation, both on 
the theoretical and the experimental side. 
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A Validation of our simulation 

In this appendix we present the validation of our simulation of the ATLAS 0-lepton and multijets 
searches. For each of the signal region in those searches we determine the 95% C.L. exclusion 
region in the tan /3 = 10, = 0, yU > slice of the CMSSM and compare to the corresponding 
exclusion region obtained by ATLAS. Such comparisons are shown in Figures A.l and |A.2 



for the 0-lepton search and the multijets search, respectively. It can be seen that, for each 
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Figure A.l: Comparison of our 95% C.L. exclusion contours with those of ATLAS in the tan /3 = 10, Aq 
^> shce of the CMSSM for the 0-lepton search. 
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Figure A. 2: Comparison of our 95% C.L. exclusion contours with those of ATLAS in the tan (3 = 10, 
fj, > slice of the CMSSM for the multijets search. 
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signal region, our exclusion contour is similar to the ATLAS one so we conclude that our 
approximation is reasonable. 
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